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ABSTRACT 

In the paper by Barta et al. (ApJ, 2010) the authors addressed some open questions of the CSHKP 
scenario of solar flares by means of high-resolution MHD simulations. They focused, in particular, 
on the problem of energy transfer from large to small scales in decaying flare current sheet (CS). 
Their calculations suggest, that magnetic flux-ropes (plasmoids) are formed in full range of scales by 
a cascade of tearing and coalescence processes. Consequently, the initially thick current layer becomes 
highly fragmented. Thus, the tearing and coalescence cascade can cause an effective energy transfer 
across the scales. In the current paper we investigate whether this mechanism actually applies in solar 
flares. We extend the MHD simulation by deriving model-specific features that can be looked for in 
observations. The results of the underlying MHD model showed that the plasmoid cascade creates a 
specific hierarchical distribution of non-ideal/ acceleration regions embedded in the CS. We therefore 
focus on the features associated with the fluxes of energetic particles, in particular on the structure 
and dynamics of emission regions in flare ribbons. We assume that the structure and dynamics of 
diffusion regions embedded in the CS imprint themselves into structure and dynamics of flare-ribbon 
kernels by means of magnetic-field mapping. Using the results of the underlying MHD simulation we 
derive the expected structure of ribbon emission and we extract selected statistical properties of the 
modelled bright kernels. Comparing the predicted emission and its properties with the observed ones 
we obtain a good agreement of the two. 

Subject headings: Sun: flares — Sun: magnetic reconnection — Sun: electron acceleration 



1. INTRODUCTION 

For many years the common picture of solar erup- 
tions and fla r es is based on the CSHKP model (see 
Ma gara et ail I1996L and references therein). It in- 
volves reconnection in a global vertical flar e current layer 
forme d behind ejected flux- rope/filament (|Lin fe Forbesl 
I2000L and references therein). The standard model is 
in good agreement with observed large-scale dynamics of 
eruptive events. Several question, however, remain open. 
In particular, it is not clear, by which mechanism the en- 
ergy accumulated in relatively large-scale (~ 1000 km) 
structures of magnetic field associated with the flare CS 
is transfered towards the small dissipation scales. In- 
deed, the magnetic diffusion in the almost collisionless 
solar coronal plasma is an ess entially kinetic process at 
scales of the order di — c/ui p i (Buchner 2006). Another 
open question is the relation between the well-organised 
dynamics of solar eruptions observed at large scales and 
the H XR and radio signatures of fragmented energy re- 
lease (jAschwandenl 120021 ; iKarlickv et all l2000f ) . And fi- 
nally, the CSHKP model has b een put in question by 
many authors (e.g. IVlahosl [2007'). in particular because 
its apparent incapability to explain the large fluxes of 
accelerated particles inferred from HXR observations. 

Recently - in order to address the issue of en- 
ergy transfe r acro ss the broad range of MHD scales - 
iBarta et al.l ([2010) [in further text referred as Paper I] 
investigated magnetic reconnection in an extended cur- 
rent layer by means of high-resolution MHD simula- 
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ti qn. Inspired by the "frac tal reconnection" conjecture 
of iShibata fe Tanumal (|2001l ) the authors of Paper I per- 
formed a 2.5D numerical simulation covering large range 
of MHD scales. They conjectured, that - in analogy with 
the vortex-tube cascade in fluid dynamics - a cascade of 
magnetic flux-tubes from large to small scale can pro- 
vide the mechanism for energy transfer across the scales. 
Two mechanisms of fragmentation were identified in Pa- 
per I: (1) The tearing c ascade - in line with conce pt of 
fractal reconnection by IShibata fc Tanumal (|2001[ ). de- 
veloped rec ently into the t h eory of ch ain plasmoid in- 
stabili ty by lLoureiro et al.1 (|2007| ) and lUzdenskv et al.1 
(l2010f). supported further by numerical MHD sim u lation s 
bv iBhattachariee et al.l (12009D. ISamtanev et all (|2009T) . 
and iHuang fc Bhattacharjed (|2010f ) and. (2) ambient- 
field driven coalescence of flux-ropes/plasmoids leading 
to formation of transversal current sheets subjected fur- 
ther to the same chain of processes of cascading fragmen- 
tation. Extrapolation of the results leads to the picture 
in which the plasmoid cascade continues down to the 
kinetic scale where actual magnetic dissipation and par- 
ticle acceleration occu rs, most likely via t he kinetic coa- 
lescen ce of plasmoids (jDrake et al.l l2005: Ka rlicky et al.1 
l2010f h Nevertheless, other processes can appear in the 
range between MHD and plasma kinetic scales, e.g., the 
Hall type of reconnection - s ee recent simulation s by 
iShepherd fc Cassakl ((20Tol ) and lHuang et all ll2010h . 

As it has been shown in the Paper I, cascading re- 
connection forms multiple non-ideal regions (at the res- 
olution limit - see discussion in Paper I) hierarchically 
distributed in the fragmenting current layer. As a re- 
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suit, the open questions of energy transfer, fragmented 
vs. organised flare energy release, and particle accel- 
eration seem to be closely related to each other in the 
presence of cascading reconnection. 

The question arises, whether this mechanism found by 
simulations is relevant for actual solar flares. We address 
this topic in our present paper. For this sake we derive 
critical signatures specific to the model developed in Pa- 
per I, that allow comparison with observations. It is cur- 
rently impossible to observe small-scale magnetic struc- 
tures, whose formation in the fragmented flare current 
layer is predicted by the cascading-reconnection model, 
directly. Almost all information on the impulsive phase 
of flares comes from the radiation emitted by acceler- 
ated particles. Therefore we concentrate on such predic- 
tions of the MHD model that are connected with specific 
distribution and dynamics of acceleration regions in the 
fragmenting flare current layer. 

In particular, by magnetic mapping of non-ideal re- 
gions embedded in the current sheet to the photo- 
sphere we derive expected consequences of the distribu- 
tion and dynamics of acceleration regions for flare rib- 
bons. Using the results of underlying MHD model pre- 
sented in Paper I we qualitatively derive the structure 
and dynamics of expected ribbon emission. We com- 
pare the modelled ribbon structure with observations di- 
rectly and by means o f its statistical properties used by 
INishizuka et ail f|2009T ). 

The paper is organised as follows: First, in Section [2TTT 
we briefly summarise the results of the underlying high- 
resolution MHD simulation of cascading reconnection. 
Then, in Section 12.21 we present the extension of un- 
derlying MHD model utilised to produce model-specific 
results/consequences in the form able to be directly com- 
pared with observations. In particular, we study the 
distribution and dynamics of dissipative/acceleration re- 
gions, represented by the embedded X-points, and their 
magnetic mapping to the chromosphere/photosphere of 
the Sun. For obtained dynamics of acceleration regions 
and their magnetic foot-points we qualitatively model the 
expected emission of flare ribbons (Section [3]). We com- 
pare the emission structure of modelled and observed rib- 
bons. Furthermore, we extract selected statistical prop- 
erties of emission kernels embedded in modelled ribbons 
and compare them with the results found by analy sis of 
actual ribbon observations (Nishizu ka et al.l j2009) . Fi- 
nally, in Section |4] we discuss our results with the in- 
tention to evaluate the relevance of cascading reconnec- 
tion/fragmentation processes presented in Paper I for ac- 
tual solar flares. 

2. MODEL 

In this section we describe the procedures allowing to 
simulate observable signatures specific to cascading re- 
connection. For clarity we first briefly summarise the 
main features of current-layer fragmentation model (Pa- 
per I). 

2.1. Basic MHD model 

In Paper I the authors studied cascading reconnection 
and energy transfer from accumulation (large) to dissi- 
pation (small) scales. They used a high-resolution MHD 
model with a non-ideal, resistive term depending on the 
current-carrier drift velocity 



vu as 

f : |ud| < v cr 

v(r,t) = ^ c iM r,i)- Vcr) . M>Vcr (1) 

with the threshold v cr set to the higher, more realistic 
value corresponding to increased resolution. The authors 
applied the model to the large-scale vertical current sheet 
in solar eruptive flares. 

The simulations in Paper I have been performed in di- 
mensionless variables: Spatial coordinates x, y, and z 
has been expressed in units of the current sheet half- 
width La at the photospheric level (z = 0) and the 
time has been normalised to the Alfven transit time 
r A = L a /Va,o, where V A ,o = B /^/fi p is the asymp- 
totic value (x — > oo, z = 0) of the Alfven speed at t = 0. 
In order to relate the model to real solar conditions ap- 
propriate scaling has been adopted in Paper I. Based on 
the gravity-introduced scale-height relation La = 600 km 
has been found. Due to this relation both the dimension- 
less and SI units can be used further in our present paper. 

In Paper I the authors described how tearing and 
driven coalescence instabilities lead to the fragmentation 
of an originally unstructured flare current layer. As a 
consequence a cascade of magnetic flux-ropes/plasmoids 
is formed from large to consecutively smaller scales. The 
structuring of the current layer leads to formation of 
multiple thin embedded current sheets. They can host 
non-ideal regions where particles can be accelerated. En- 
hanced numerical resolution revealed the structure of 
dissipation regions: They form many thin channels of 
non-zero magnetic diffusivity. The spatial distribution 
and dynamics of these diffusivity channels has been fur- 
ther studied via tracking of the X-points associated with 
the non-ideal regions. Analysis shows hierarchically- 
structured grouping of the dissipation channels and in- 
tcrmittency in their life-times and times for which they 
are magnetically connected to the photosphere. 

The structure, distribution and dynamics of the non- 
ideal regions should be reflected in specific observable 
features that we are going to derive in the following sub- 
section. 

2.2. Derivation of model- specific observable features 

Unfortunately, it is impossible to directly measure the 
magnetic field or current density in the coronal current 
sheets (and at all not with the resolution reached in the 
simulation done in Paper I). Hence, we need some in- 
direct, more subtle comparison between observed and 
modelled quantities. A possible way which may lead to 
indirect indication of curr ent-sheet fragmentatio n has re- 
cently been presented by INishizuka et al.l (|2009f ) . These 
authors study the structure of emission in flare rib- 
bons, namely the distribution and dynamics of embedded 
bright kernels. They conjecture possible relation of the 
resulting power-law distributions found in their statisti- 
cal analysis of ribbon- k ernel properties with the c oncept 
of fractal current sheet (jShibata &: Ta numa 200ll) . 

We follow this idea from the other end. We start 
with the structure, distribution, and dynamics of dif- 
fusion/acceleration regions that are specific to the 
cascading-reconnection model as they has been found in 
Paper I. Then, we relate these features to the structure 
and dynamics of emission in flare ribbons. We proceed as 
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follows: As shown in Paper I (Fig. 6), the reconnection 
X-points are associated with thin channels of magnetic 
diffusivity. Hence, we take the X-points as geometric rep- 
resentatives of dissipation regions. In these regions elec- 
trons can be accelerated, e.g., by the DC electric fields. 
However, also diffe rent acceleration m echanisms has been 
proposed - e.g. iDrake et al.l (|2005l) suggest Fermi-type 
acceleration in contracting plasmoids. We can include 
this type of acceleration into account by following con- 
sideration: Extrapolating the results of numerical sim- 
ulations of reconnection cascade to smaller scales one 
can imagine, that each dissipative region surrounding 
studied X-point in fact contains many unresolved small- 
scale magnetic islands. Electron acceleration in these 
no n-ideal regions is t hen performed (in line with ideas 
bv IDrake et all I2005T) via coalescence and shrinkage of 
these micro-plasmoids. Anyway, electrons accelerated by 
Fermi-type mechanism inside the plasmoids are trapped 
and can be released to the open field (which connects CS 
to chromosphere) again in the vicinity of the X-points, 
where they become demagnetized. For further drelated 
discussi on on elect r on ac celeration in magnetic islands 
see, e.g. lOka et al.1 (goTS). 

In our study we focus on electrons as we are interested 
in optical and UV/EUV chromospheric response ascribed 
to electron beams. These are then transported along the 
magnetic field lines either up-ward to the solar corona 
or downward, until they reach dense layers of the solar 
atmosphere, where their energy is thermalized. Here we 
concentrate on the down-ward transported electrons. In 
order to study the positions and dynamics of the points, 
where these electrons reach the chromosphere, i.e. the 
expected bright kernels, we map found X-points to the 
bottom boundary of the simulation box using magnetic 
field. The magnetic field lines that go through given X- 
point represent magnetic separatrices in our 2D model. 
Therefore, we are seeking for intersections of magnetic 
separatrices with the bottom boundary. In the following 
we will refer to this footpoints at the bottom boundary 
as to kernels. Thus, from known position of all X-points 
(Fig. 7 in Paper I) and knowing the magnetic field we 
can calculate the positions Xk{t) of all kernels k for each 
recorded time step t of the simulation. Note, that due 
to the 2D geometry of the model the foot-point/kernel 
position is uniquely given by its x-coordinate (see Paper I 
for details). 

The calculations of the chromospheric emission in a 
certain spectral line (e.g. H Q or C iv) in response to 
the bombardment of electrons acceler ated in a distant 
reconnection region is a difficult task (Kasp arova et all 
120091: iVaradv et al.ll2010() . In order to obtain an at least 
qualitative output that could be later compared with ob- 
servations we use for ribbon-intensity distribution the ex- 
pression 

T i .n j ( ( x ~ x k(t)f\ 
I k (x, y, t) = I exp I — I 

I(x,y,t) = J2l k (x,y,t) , (2) 

k 

instead of description of all the complicated processes of 
electron transport, energy deposition, chromospheric re- 
sponse, and radiative transfer. The intensity I is summed 
over all kernels k; y is the second (however invariant in 



our model) coordinate in the photosphere. The kernel 
size is set to A = O.ILa- Since the model is capable nei- 
ther to estimate the electron-beam flux nor the energy 
thermalized in the kernel we set the intensity scale to 
Iq = 1. We will return to this point later in the discus- 
sion. 

Let us now apply this procedures to the results of sim- 
ulations presented in Paper I. 

3. RELATION TO OBSERVATIONS 

Fig. [1] depicts the above-described procedure of map- 
ping of X-points to the photosphere along magnetic sepa- 
ratrices and flare-ribbon emission calculation. Obtained 
results are shown for t = 391. 5ta- Panels (a) - (c) show 
how three dissipative regions/X-points embedded in the 
fragmented current layer are magnetically connected to 
the bottom boundary and mapped by separatrices to the 
positions x f=s ±2800 km (one kernel) and x ~ ±3750 km 
(two mutually close kernels). 

A view at the flare ribbons obtained by applying 
the simple model described by Eq. ^ is presented in 
Fig. [ljd). As the figure indicates, in some cases one 
can observe even multiple pairs of flare ribbons. This is 
clo se to the idea o f chrom ospheric re-brightening studied 
bv lMiklenic et~"aTl (|2010h . The inner separatrix associ- 
ated with the internal pair of ribbons is connected to the 
X-point which appeared temporarily in the transversal 
current sheet formed between merging plasmoid and the 
loop arcade (see also Fig. 4 in Paper I). Figs.QJe) and (f) 
show the emission intensity profile along the a;-axis cal- 
culated according to the model in Eq. |2]). The enlarged 
detail shown in Fig. QJf) reveals an internal double-peak 
structure of the outer ribbons. Note, that a similar struc- 
turing is seen in observed H a ribbons (Fig. [2]). 

In order to study the dynamics of kernels and associ- 
ated X-points during the entire recorded interval of evo- 
lution (300ta — 400ta) we tracked the moving positions 
of all X-points and the magnetically-associated (mapped) 
kernels. Kinematics of the kernels (foot-points of separa- 
trices) is depicted in Fig [3] Panel [Sja) presents a global 
picture of kernel dynamics over the entire recorded inter- 
val. Because of the CS symmetry only the right ribbon 
positions are displayed. The spatial coordinate x, the 
distance from the polarity- inversion line (PIL), is limited 
to the interval where the kernels actually occur. Fig[3](b) 
shows a detailed view of the area selected in panel (a). In 
order to reach higher zoom the global trend (motion) of 
the group of kernels (i.e. the increasing ribbon separation 
as the flare proceeds) has been subtracted by applying 
the transformation x/La — x/La — 0.0125(£/ta — 300). 
X-points that are separated only by tiny magnetic is- 
lands can map practically to the same kernel. They 
then cannot be distinguished even in the zoomed dis- 
play. Therefore, for each time we take a set of kernel 
positions (its size varies with time) and display every 
kernel in different colour. Fig. [3] shows various life-times 
of the kernels and other aspects of their dynamics, like 
bifurcation/merging. Note, that the life-time of the ker- 
nel is not determined solely by the life-time of the asso- 
ciated distant X-point but also by the processes in the 
current layer that can change the magnetic connectivity 
of diffusion regions embedded in the current layer to the 
photosphere (Fig. 7 in Paper I). 

We have further extracted selected characteristics of 
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Fig. 1. — Model-specific consequences of cascading magnetic reconnection. (a) Global 3D magnetic and current-density structures at 
t = 391. 5ta- The magnetic field (black lines), the separatrices that map the diffusive regions to the bottom boundary (green lines), and the 
current density (red colour scale) are shown, (b) Projection of (a) int o the xz-pl&ne. The rr-axis shows positions both in the units of La 
(top) and in kilometres according to the scaling adopted in Section [Zl] (c) Enlarged view of the selected rectangle reveals double-structure 
of the outer pair of separatrices hitting the bottom boundary near x ± 3750 km. (d) Modelled view at two pairs of flare emission ribbons 
(inverted colour scale, darker means higher intensity here), (e) Modelled emission profiles across the ribbons, (f) Detailed view of the 
outer-ribbon profile reveals its internal double-peaked structure. 



the modelled kern el dynamics and comp are them with 
the observations of iNishizuka et all (|2009) . For this sake 
we processed our model results in the same manner as 
those authors by binning photosphere along the x-axis. 
For the bin size we chose Ax = 0.05La- For each time 
we integrated the intensity given by Eq. §2§ over each 
bin. Thus we obtained a proxy for "light-curves" of all 
bins. Fig. |4] shows three examples of such light-curves 
for bins centred around positions x = 4.58La (2750 km), 
x = 4.93L A (2950 km), and x = 5.73L A (3435 km). 
Narrow peaks coming from the short-living X-points are 
superimpos ed over the longer liv ing bright kernels (c.f. 
Fig. 2( b) inlNishizuka et alj|2009D . 

As in INishizuka et al.l (|2009l ) , we then determined the 
peaks at each light-curve and recorded their time of oc- 
currence and their maximum intensity. We calculated 
the statistical distributions of the peak intensities and 
the time-intervals between two subsequent peaks. The 
results of this analysis are presented in Fig. [5] Despite 
of the poor statistics (for our ID mesh we identified only 
68 peaks) we obtained spectral indices s — —1.48 for the 
peak-intensity (panel E)(a)) and s = —1.73 for the time 
interval between subsequent peaks (Fig. [5jb)). 

4. DISCUSSION AND CONCLUSIONS 

The 'standard' CSHKP picture of solar eruptive flares 
involves the formation of a global flare current sheet (CS) 
behind a ejected filament followed by magnetic reconnec- 
tion in this CS. Nevertheless, deeper study of this sce- 
nario invokes further questions: 1) How is the energy 
accumulated in flares in relatively large-scale structures 
(global flare current layer) transferred to the dissipative 
(kinetic) scales, 2) How can we observe regular dynam- 
ics in eruption/flare on large scales, just in line with the 
CSHKP scenario, and the signatures of fragmented en- 
ergy release at the same moment, and 3) How can the 
large fluxes of accelerated electrons inferred from HXR 
observations be reconciled with a single, relatively small 
diffusion region assumed in classical CSHKP model. 



The simulations presented in Paper I showed, that 
a cascading reconnection resulting from the formation 
and interaction of plasmoids/flux-ropes can address these 
three fundamental questions at once and that it might 
provide viable scenario for magnetic energy dissipation 
in large-scale systems, like solar (eruptive) flares. 

In the present paper we evaluated the relevance of 
this model for actual solar flares. For this sake we de- 
rived observable model predictions that are specific for 
the cascading-reconnection scenario, and searched for 
the predicted features in observed data. Since multi- 
ple acceleration regions embedded in the global flare 
current layer are - via magnetic field-line mapping - 
inherently connected with the structured emission in 
the flare ribbons, detailed study of the ribbon kernels 
might bring some evidence for the cascading processes 
in the current layer above the flaring region. This 
idea is illustrated by Fig. [1] which shows the struc- 
turing of modelled emission in ribbons. Fig. [1] fur- 
ther indicates, that the larger-scale plasmoid interact- 
ing with the flare-loop arcade can temporarily form a 
second pair of ribbons. The plasmoid/loop-arcade in- 
teraction may play significant role in flare dynamics. 
Another indirect evidence of such interaction based on 
analysis of series of X-ray images of the limb flar e has 
been presented bv iKolomanski fc Karlickyl (|2007| ) and 
iMilliean et all (I2010D . 

Similar structuring of H a emission as presented in 
Figs.QTd)-(f) ig 5 indeed, visible in Fig. [2] Yet another in- 
teresting feature in Fig. [2] is the dark filament-like struc- 
ture (indicated by an arrow). Since the original large- 
scale filament whose eruption initiated the flare was al- 
ready far away at the time of observation, the absorp- 
tion feature can be interpreted as a manifestation of one 
of the secondary plasmoids formed in the global current 
layer in line with the presented cascading-fragmentation 
scenario. The secondary plasmoids represent enhanced- 
density structures. Density increase in connection with 
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Pixel position along the path 



Fig. 2. — Observation of flare H a emission ribbons (a) and 
profiles of H a intensity (b) along selected paths (red lines in the 
panel (a)). The ribbon emission is structured (double-peaked?) 
— c.f. Fig. [T] The green arrow points to the filament-like dark 
structure which projects to the space between flare ribbons. 



consequent faster radiative cooling might lead to the de- 
tectable H Q absorption. The secondary plasmoid can be 
consequently manifested as a darker feature at the back- 
ground of relatively brighter chromosphere. 

A sophisticated stud y of the structu r e of flare rib- 
bons was presented by iNishizuka et al.l (|2009f) . These 
authors studied the statistical properties of emission ker- 
nels found inside the ribbons and they found a power- 
law distribution of the studied kernel characteristics. 
A power-law distribution ca n be a signature of self- 
organ ised criticality evolution (|Aschwandenll2002l IVlahosI 
2007 f). It can also indicat e a fractal current sheet decay 
(Shib ata fc Tanumall200il ). The latter scenario is in fact 
favour ed by the subsequent analysis of INishizuka et all 
(2010) which relates the energy release measured by 
means of HXR flux and ejection of multiple plasmoids 
formed in the fragmenting current sheet. 
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x/L A -0.0125*(t/r A -300) 

Fig. 3. — Kinematics of flare ribbon kernels obtained by magnetic 
mapping of dissipation regions in cascading reconnection down to 
the photosphere, (a) View of all kernel positions in t=300 ta — 
400 t\. Abscissa positions in kilometres (upper scale) and in units 
of La- (b) Detailed view of the area selected in (a) (the blue 
rhomboid). The global motion of the group of kernels has been 
subtracted in order to reach larger zoom. For each time instant 
each kernel is depicted by unique color in order to distinguish be- 
tween close (almost overlapping) kernels/foot-poits. 



In our present paper we have established the connec- 
tion between the hierarchical distribution of diffusion re- 
gions formed by cascading reconnection and the statis- 
tical properties of emission kernels from first principles. 
We have studied the dynamics (positions, and the mech- 
anisms and times of creation/annihilation) of O- and X- 
type points in the current layer as well as of the kernel- 
points associated with the diffusion regions (X-points) 
by means of magnetic field-line mapping. We then asso- 
ciated each found kernel with the emission according to 
Eq. ((2|). We chose this simple model as we do not have 
any information about the accelerated electron fluxes ab- 
sorbed in the kernel, nor can we easily calculate the chro- 
mospheric response in terms of (e.g., H Q ) emission. Our 
choice of the unified peak intensity Iq = 1 for all kernels 
can be justified following way: After all, the particle ac- 
celeration is performed by kinetic processes at the small- 
scale end of the c ascade. Suggested mechanism (see e.g. 
iDrake et al.ll2005l) involves, e.g. kinetic coalescence and 
shrinkage of small-scale magnetic islands/plasmoids. In 
the dissipative range of scales the slope of the energy 
spectrum (see Fig. 1 in Paper I) is very steep. Hence, the 
size of plasmoids at kinetic scale can be expected more 
or less uniform as indicated also by PIC simulations. In 
such case the elementary injections of accelerated parti- 
cles are, perhaps, roughly of the same magnitude. The 
intermittent, spiky structure of observed emission pro- 
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Fig. 4. — Modelled light-curves for three selected bins at positions 
x = 4.58L A (2750 km), x = 4.93L A (2950 km), and x = 5.73L A 
(3435 km). 



files is then the result of o verlapping of many of these 
elementary injection peaks (|Aschwandenll2002l) . Indeed, 
we observe this feature in examples of the modelled light- 
curves obtained at fixed points at the bottom boundary 
(Fig. |4|). Despite of the intensity produced by the ele- 
mentary particle injection being normalised to Jo = 1, 
the overlapping elementary peaks produce spiky inten- 
sity profiles ranging from / = to I > 3. 

In or der to directly compare our results with the find- 
ings of iNishizuka et al.l (|2009f ) we have performed the 
same statistical analysis of the peaks in our artificial 
light-curves. We have found the statistical distributions 
of the emission peak intensities and of the time-lag be- 
tween subsequent emission peaks (Fig. [5]). Because of 
the 2-D geometry of underlying MHD model we have 
obtained 1-D distribution of emission kernels. Conse- 
quently we have found smaller number of the peaks. De- 
spite of the resulting poor statistics we could conclude, 
that the slopes of the power-law distributions obtained 
by our model agre e surp risingly well with those found 
bv INishizuka et all ()2009D for observed light -curves: For 
the statistical distribution of the peak intensities we have 
found spectral index s = -1.48 ()Nishizuka et al.l (|2009l ) 
obtained s — —1.5) and for the time-la g between peaks 
we ob tained s — —1.74 (s = —1.8 in Nishizuka et al.l 
(2009)). We are aware that our statistics is, however, 
rather poor as it operates with only 68 peaks found in to- 
tal. Thus, it would be premature to take this agreement 
as conclusive. Nevertheless the result indicates that the 
statistical properties of emission kernels resulting from 
cascading reconnection and of those observed are - at 
least - not in direct contradiction. 

Results of presented comparison between the modelled 
features attributed specifically to cascading reconnection 
and their observed counterparts supports - in our view - 
the relevance of cascading reconnection processes studied 
in Paper I for reconnection in solar atmosphere. Agree- 
ments found between the statistical properties of mod- 
elled and observed emission kernels invokes, however, fur- 
ther questions. The underlying 2.5D MHD model ig- 



a) 



-1.48 



10.0 



1 .0 r 




-0.4 -0.2 0.0 0.2 0.4 0.6 

log(l) 



100.0 



b) 



10.0 



1 .0 



0.1 



-(At) 



1.74 



-0.4 -0.2 



0.0 0.2 

log(At/r A ) 



0.4 



0.6 



Fig. 5. — Statistical distribution of the modelled intensity (a) 
and the time-interval between subsequent emission peaks (b) for 
accumulated light-curves from all mesh-boxes. N indicates the 
number of peaks with the property within the bin. 



nores structuring of the current sheet in the y-direction 
(along the PIL). This result might indicate that modes 
with k y > (in the off-plane direction) may not play 
too significant role in current-layer fragmentation in so- 
lar flares, or at least, that they do not reach the smallest 
scales to substantially influence the dynamics (life-time) 
of kinetic-scale diffusion regions. This issue can be ad- 
dressed, however, only within the framework of a full 3D 
MHD model. 
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